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1. Introduction 

In the papers [3] and [5] the Schur transformation for gcnerahzed Nevanhnna func- 
tions with a reference point zi in the open upper half plane was considered. An 
analogous transformation for Nevanlinna functions (for the definition of a Nevan- 
hnna function see Section [2]) and for the reference point oo is defined in jl] Lemma 
3.3.6], see [3]. This transformation or a simple modification of it we call here the 
Schur transformation for Nevanlinna functions, and it is the starting point for the 
present paper. To give more details, we consider a Nevanlinna function n which 
has for some integer p > 1 an asymptotic expansion of order 2p + I at oo, for 
example 

f ^ So Si S2p f I \ . , . 

"(^) = -7-^ +''[^2^)^ z = 12/, y ^ ±oo. (1.1) 

The Schur transform n of n is the function 

n(.):=-^-z+fi; (1.2) 
n[z) So 

the relation between n and n can also be written as 

= si • 

z 1- niz) 

So 

The transformed function n =: ni is again a Nevanlinna function, but in general 
with an asymptotic expansion of the form (|l.ip of lower order 2p — 1, and if p > 1 
the Schur transformation can be again applied to ni etc. As a result we obtain a 
finite sequence of Nevanlinna functions rii ~ n,n2 = rii, . . . ,np = n^^i; this is the 
sequence of functions that appears in the asymptotic expansion of n by continued 
fractions, see [U Section 3.3.6]. 

The transformation p.2p is closely related to the finite Hamburger moment 
problem. We recall that the Nevanlinna function n with an asymptotic expansion 
p.ip admits an integral representation 

J — oo ^ ^ 

where cr is a bounded non-decreasing function on R. The coefficients Sj in (jl.ip 
are the moments of the function a: 

t^da{t), j = 0,l,...,2p. (1.3) 

The moment problem we have in mind is the problem to determine all Nevanlinna 
functions n with an expansion (|l.ip and given coefficients s_,, j = 0, 1, . . . ,p, see 
[H Theorem 3.2.1]. 

An essential feature in our studies are operator representations or so-called 
realizations of Nevanlinna functions, see [T3j, [THji and [5]. In fact, if the Nevanlinna 
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function n admits an asymptotic expansion (jl.ip its operator representation takes 
the simple form 

n{z) ^ {{A- z)-^u,u), z^z*, (1.4) 

with some Hilbert space Ti. with inner product (■,■), u € Ti., and a self-adjoint 
operator A in H. We study the corresponding operator representation of the Schur 
transform n, and also of the functions 712, . . . ,7ip. For example, the function n 
admits an operator representation of the form (|1.4p with a Hilbert space Ti, an 
operator A, and an element u which are the orthogonal complement of the element 
u in Ti., the compression of A to Ti., and a multiple of the projection of Au onto 
Ti., respectively. After applying the Schur transformation p times, the resulting 
function Up admits an operator representation of the form (|1.4[) with the space 

Hp =HeHp, Hp -.^ span{u, Au, . . . , AP~^u}, 

the operator that is the compression of A to this space, and an element which is 
a multiple of the projection of A^u onto Hp. 

Since Hp is obtained by subsequent application of fractional linear transfor- 
mations of the form (|1.2p . there is a fractional linear relation between the function 
n and the transformed function Up. We derive an explicit form for the defining 
2x2 matrix function V of this relation in three ways: By calculating the resolvent 
of the operator A in its 2 x 2 block matrix operator form corresponding to the 
decomposition H = Hp® Hp, by means of the description of all generalized resol- 
vents of a certain symmetric restriction of A with defect one in the space Hp+i, 
and via reproducing kernel methods using the non-negative Nevanlinna kernel 

L^{z,w) = z,w^£\R, z^w*. 

z — w* 

For the Nevanlinna function n with an asymptotic expansion (|l.ip . polyno- 
mials Cj and dj, j = 1,2, ... ,p, of first and second kind can be defined by the 
well-known formulas, see [U Chapter I] . Recall that Cj is a polynomial of degree 
j, and that dj is a polynomial of degree j — 1. We show that the polynomials ej 
of first kind of the transformed function n coincide, up to constant factor, with 
the polynomials dj^i{z) of second kind for the given function n, whereas the poly- 
nomials dj of second kind for n are linear combinations of e^+i and dj+i. As a 
consequence, the polynomials of second kind for n are orthogonal with respect to 
the measure generated by the non-decreasing function a in the representation of 
the form (|1.4p of the Nevanlinna function n; in this statement n can be replaced 
by the function ~l/n. As in the classical moment problem, the 2x2 matrix func- 
tion V, which determines the fractional linear relation between n and rip, can be 
represented by the polynomials of first and second kind. 

A short synopsis is as follows. The Schur transformation is defined in the 
next section. We start with weaker forms of the asymptotic expansion (jl.ip , for 
example 

f \ So Sl / 1 \ . 

n(z) = 2+°(~)' 2 = iy'y^±oo; 
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and consider also a weaker form of the Schur transformation. In Section[3]we men- 
tion three concrete forms of the operator representation of n. The basic resuh of 
this section is Theorem l3. II which describes the operator model for the transformed 
function. Higher order approximations and the corresponding polynomials of first 
and second kind are introduced in Section 31 In the operator model an asymptotic 
expansion can be characterized by the fact that u G domAP. The main result 
of this section is the relation between the polynomials of first and second kind 
of n and n which was mentioned above. The reduction via a p-dimensional sub- 
space, that corresponds to p subsequent applications of the Schur transformation, 
is given in Section [5] by means of a block operator matrix representation of A. In 
Section [6] the corresponding transformation matrix V is expressed in terms of the 
polynomials of first and second kind. Although the final formulas are well known 
(see for example [T]) this approach seems to be new. 

In Section [71 applying the theory of w-resolvcnt matrices, we derive a rep- 
resentation of a transformation matrix in an explicit form by means of the given 
moments; it corresponds to Potapov's formula for the solution matrix of the Nevan- 
linna - Pick problem, compare also [2| . Finally, in Section [8] we explain the con- 
nection between n and Up through some basic results from the theory of resolvent 
invariant reproducing kernel spaces, and give another proof for the representation 
of the transformation matrix by orthogonal polynomials. 



2. The Schur transformation 

1. A Nevanlinna function is a complex function n which is defined and analytic in 
the upper half plane C+ and has the property 

2 e C+ =^ Im n{z) > 0. 

We always suppose that n is extended to the lower half plane C~ by the relation 

n{z)^n{z*)*, zeC-, (2.1) 

and to those points of the real axis into which it can be continued analytically. 
The set of all Nevanlinna functions is denoted by Nq. Recall that n G Nq if and 
only if n is analytic in C \ M and the kernel 

£„(^,^;)= Z,weC\R,Z^W% 

z ~ w* 

is positive definite. 

Let n G No and consider the following properties of n: 

(lo) = + 

(2o) n(2) = -^ + 0(^l), 
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(3o) n{z) = -— -^+o( — 

where here and in the foUowing, the hmit relations are understood to hold for 
z ±00 along the imaginary axis. The assumption (j2.ip implies that sq and si 
are real numbers. Evidently, (3o) =^ (2o) =^ (lo). The function n satisfies the 
assumption (lo) if and only if it belongs to the class (Ra) of [12], which means 
that it admits an integral representation 

.+00 ^ 

n{z) = / da{t), zeC\R, (2.2) 

J —00 t ^ z 

where cr is a bounded non-decreasing function on R. Then 

da{t) = So, 

hence sq > 0, and if so = then n{z) = 0. With the representation (|2.2[) of n the 
assumption 



\t\da{t)< 00 (2.3) 
implies that (3o) is satisfied. Indeed, (|2.3p implies that 

/oc 
t da{t) 
-00 

exists and with z ~ iy 



V Z Z^J j_^t-Z J-aot +V 

The assumptions (Iq), (2o), and (3o) are all different. To see that (Iq) 7^ (2o) 
we show that if n g Nq has the representation (j2.2p with supper C [0,oo) and 

/>oo />oo 

/ da{t) < 00, / tda{t) = 00, 
JO Jo 

then (2o) does not hold: Let c > be given arbitrarily (large) and choose K > 

such that / td(j(t) > c.li y is chosen large enough then for < t < X we have 
Jo 



and hence 



t2 + y2 - 2' 
^ ■ da{t) > 



,2 



^2 + y2 W - 2 ' 



and therefore, with z ~ iy, 

t ^ , , f°° -yH + iyt^ 
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which imphes that (2o) does not hold. Thus, for example, the function 

-1 1 dt 



n{z) 



z — \/—z 



t- Z TT{t+l)Vi 



satisfies (lo) but not (2o). 

Let n be the Nevanlinna function, defined in the upper half plane by 

z + -f + f{z) 

where sq is a positive real number, 7 is a complex number with Im7 > 0, and / 
is a Nevanlinna function such that f{z) ~ o(l). It has the properties 

lim 2^ (iT-i^) + ) = 7S0, hm ('^(^) H — ~ ) — 7*^0, 

2=iy,y— >oo V z / z=iy,y—> — oc \ z ) 

and hence n satisfies (2o) but, since the two limits arc different (and non-real), it 
docs not satisfy (3o). 

Instead of the assumption (3o) also the assumption 

"«-?-i^ + o(i) (^.^) 

seems reasonable. However, according to [13l Bemerkung 1.11], (|2.4p implies the 
existence of a real number si such that 

n(.) = -^-^-^+o(l); (2.5) 
z z-^ \z J 

this relation will be considered in Section [3] as assumption (li). The implication 
(|2.4p =^ (|2.5p can also be seen from the integral representation (|2.2|) of n: (|2.4|) 
implies 

1 +irf,(i) + fi\ 0(1), 

\. ^. t - z z z^ 



and hence with z = iy 

yH{t + iy) 



, t2 + y2 

Taking the imaginary part we see that 



da{t) + iysi ~ 0(1), y ^ co. 



yH 



and taking the real part we see that there exist real numbers C and j/o such that 



■da{t)<C, y>yo. 



This implies that 



/oo 
da{t) < 00, 
-00 
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hence 

/oo 
\t\ da{t) < 
'OO 

and 



/OQ 
t da{t). 
-oo 



Now (|2.5[) easily follows from the integral representations of n and the expressions 
for the real numbers sq, si, and S2- With z = iy we have 



z 



-.f , . So si S2 \ r , , , f°° it^y , , , 



+ 

2. Now we define the basic transformations considered this paper. 

Definition 2.1. If n e No satisfies the assumption (Iq) or (2o), the Schur type 
transform n of n is the function 

n(z) = ^-z, (2.6) 
if n e No satisfies the assumption (3o) the Schur transform n of n is the function 

ni,) = ^ (2.7) 
n(z) So 

The difference between the formulas (|2.6p and (|2.7p is just in the additive 
real constant si/sg: under the stronger assumption (3o) this constant assures that 
the transform tends to zero if z tends to ±oo along the imaginary axis, see (j2.10p 
below. 

The relations (j2.6p and (|2.7p can also be written as a first step of a continued 

fraction expansion 

/ N So , , So 

n[z) = ^, , , or n[z) = ^ , , . 

^ ' z + n{z) ^ ' z- ^+n{z) 

Theorem 2.2. The following equivalences hold: 

n e No and satisfies (lo) n e No, n(z) = o(z), (2-8) 

n e No and satisfies (2o) n e No, n{z) = 0(1), (2.9) 

n e No and satisfies (3o) n e No, n(z) = o(l). (2-10) 

Proof. We have 

~, . ^, X si So + zn{z) 

n{z) = n{z) = — . 

So n(z) 

A straightforward calculation yields 

Imz / Imn(z 



Imn(2;) = Imn(z) = . ( so "7'"^"^ - \n{z)\'^ 

|n(z)|^ \ imz 



and the estimate 



Kz)|^ = 



+ 00 



da(t) 



-oo ^ Z 



2 



dcj{t) r°° , . . _ Imn(z) 
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implies n,n£ Nq. The asymptotic properties of n follow in case (|2.8p from the 
relation 

n{z) So 



z zn{z) 
from the relation 



_ So + zn{z) 

n{z) = -- 



1=0(1), 



zO[ - 



n{z) zn{z) ' 

and for n(z) in case (|2.10|) in a similar way or from [1] Lemma 3.3.6]. 
Conversely, starting from n{z) as in (|2.8p . the relation 

n{z) 

z (n{z) + — ) = So 

V z J n{z) 

z 

implies that from n{z) = o(z) it follows that n satisfies (lo)- The corresponding 
proofs for (|2.9|) and (j2.10p are similar. □ 



3. Self-adjoint operator representations 

A function n € No admits a self-adjoint operator representation or realization with 
a self-adjoint relation A in some Hilbcrt space TC of the form 

n{z) = nizoT + {z- z*^) ((/ + {z ~ zo){A - z)-') v, v) (3.1) 

with zq an arbitrary non-real number zq and an clement u G 7i, see |13| . |10| . and 
[5]. If w is chosen to be a generating element for A, which means that 

n = span {v, {A - z)-^v | z e C \ M} 

and which is always possible, then the operator representation (|3.2p is called min- 
imal and then it is unique up to unitary equivalence. We have the following equiv- 
alences, see [15] : 

A is an operator <=^ n{z) = o(z), 

V € dom74 <=^ limj,_,oo2/Inin(i?;) < oo; 

for n € No the latter limit always exists: it is either a non-negative number or oo. 

If the Ncvanlinna function n satisfies the assumption (lo) (or any of the 
assumptions (2o), (3o)) the representation (j3.ip can be simplified to 

n{z) {{A- z)-\,u), zeC\R, (3.2) 

where A is a self-adjoint operator in some Hilbert space 7i, u € 7i, {u, u) = sq. If 
u is chosen to be a generating element for A, or equivalently, 

n = span {{A- z)~\ \zeC\R} 



9 



which is always possible, then the operator representation (|3.2p is also called min- 
imal and then it is unique up to unitary equivalence. The representation (j3.2p 
follows from (j3.ip and the above mentioned equivalences by taking u = c{A — zq)v 
with some unimodular complex number c. 

Here are three examples for a more concrete choice of the triplet Ti., A, u in 
p.2[) for the given function n E No with integral representation (|2.2[) . 

(1) Ti. = L^(a), A is the operator of multiplication with the independent variable, 
and u{t) = 1, f e R. 

(2) Ti is the completion of the linear span of the functions r^, z e C \ M: 

t — z 

with inner product defined by 

A is the operator of multiplication by t, and u{t) = 1, t e M. 

(3) Ti is the reproducing kernel Hilbcrt space C{n) with reproducing kernel 

L^{z,w) = ^(^)-^H\ z,weC\R,z^w% 
z — w* 

A is the self-adjoint operator whose resolvent [A — z)~^ is the difference- 
quotient operator i?^: 

and take u — n; this function belongs to the space C{n), since n satisfies 
the condition (Iq)- Recall that the reproducing property of the kernel i„ is 
reflected in the inner product of the space C{n): 

(/, Ln{ ■ , z)c(n) = fiz), f e C{n), zeC\R. 
That (IMl) holds follows from 

(i?,n)(C)=i„(C,^*) 
and the reproducing property of the kernel i„: 

{{A - zy^U, u)c(n) = {RzTl, n)£(„) 

= (ri, i„( • , z*))£(„) = n{z*)* = n(z). 

The unitary equivalence of the representations in (1) and (2) follows easily 
from the relation 
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and the fact that the functions r^, z ^ z* , form a total set m LF'{a). The unitary 
equivalence between the two representations of n in (2) and (3) is given by the 
mapping U : 

[/(r,) = L„(.,z*); 

in particular, we have Uu — n where u is the function u{t) = 1, < e K. The space 
L'^{cr) (or the equivalent space in (2)) wc denote also by 'H{n). We mention that the 
definition of the spaces in (2) and (3) can also be used for generalized Ncvanlinna 
functions, whereas in this case the space i^(o') need not be defined. In Sections 
|4]-[7]we will prove theorems using the representation of n in (1), in Section [8] we 
reprove some of these results using the representation in the reproducing kernel 
Hilbert space C{n). 

Since, according to Theorem 12.21 the hmctions n and n in Definition 12.11 
belong to the class No and are o(z) for z = iy, y oo, they admit again an 
operator representation of the form (|3.1[) . for example. 



n{z) = n(zo)* + {z- z*) ((1 ~ zo){A - z)-'d, (3.3) 

with a self-adjoint operator A in some Hilbert space 7i, zq an arbitrary non-real 
number, and an element v € Ti.. Clearly, as the difference between the functions n 
and n is just an additive real constant, the operator representation for n can be 
chosen the same, that is, in (j3.3p n can be replaced by n. 

Theorem 3.1. Let n e Nq satisfying the condition (Iq) and with operator repre- 
sentation p.2[l be given, and let 

. -So , Si 

n{z) = —-^ -z + — 
n[z) so 

be the Schur transform of n from (|2.7p . Then in the operator representation p.3p 



ofn we can choose "H = {w}^, A inTi as the compression of A to "H : A ^ PAU 



where P is the orthogonal projection in Ti onto TL, and the element v as 



u 



{(A-Zq) ^u,u) 
Ifn also satisfies the condition (lo) 0? then 



PiA-zo)-'u 



PAu 



Remark 3.2. The resolvent of A is given by 



This is the case when n satisfies condition (li) defined in Section |4l see Lemma |4. II 
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and 

u 

\\u\\{A- zo)-^u- {{A - za)-^u,u) — 

\\u\\_ 

\{{A-zo)-^u,u)\ 

Note that ((^1 — zo)^^u, uj = n{zo) ^ 0, otherwise n{z) = 0. 



Proof of Theorem \3A[ (1) Suppose that n satisfies (lo). Then we have = ^/sq 
and 

n{z) - n(zo)* = - Vt - ' + W + (3-4) 

where we have put r{z) := {{A — z)~^u, u). It remains to show that the expression 
on the right hand side of p.4p equals 

(z-z*) [iA-zo)iA- z)-')d,d) = {z~z*)\\d\\' + iz-z*){z~zo) ((1- z)-'d,d) . 
This is a straightforward calculation, we only indicate some formulas: 

iz-zoKA^z)-'d = J^^^f'i^^A-zr'u-iA-zor'u 

\r{za)\ \ r{z) j 

and 

^11 1 {z-zl)\\{A-z^)-^uf 

""II \r{z^Y t{z) |r(zo)P 

(2) Now assume that n satisfies (lo)- Then v G doniA C dom A and the equality 

'kll , , 1 u 



-{A- z^Y^u - — 



r(zo) ||w|| 

shows that also u S dom A. If we take m = —i^A — zq)v, then (see after (|3.2p ) n has 
the asserted representation. It remains to show that u = PAu/\\u\\. We have 

u = -{A - zq)v = -^{A - zo)u - 4^"- 
\\u\\ r{zo) 

Taking the inner product of both sides with u and using (m, u) ~ 0, we see that 

_ {{A- zq)u,u) 
r{zo) ||w||3 

and hence 

I , ((A-zq)u,u) 1 f , (Au.u) \ 1 

u= (A-z^Y- "-^ " > u=—(Au-^-^^u]=—PAu. 

□ 
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4. Higher order asymptotics. Orthogonal polynomials 

1. For n e No and some integer p> \ we introduce the assumptions 

(Ip) "(^) = "7-^ ^2^+°^;2^ 

(2p) n(z)----- ^^ + 

(3p) n{z) = - — + o ^ 

Again, (3p) =^ (2p) =^ (Ip), and by [l3l Satz 1.10] for the operator represen- 
tation the assumption (Ip) is equivalent to u G domA^. That is, for the above 
representation with the space 7i(n) the functions 



22p+2 
1 



tfc(^):=^^ A: = 0,1,, 



belong to Ti{n) and the first p of these elements, to, ti, . . . , tp_i, belong to dom A. 
Moreover, the formal relation 

2P AJ A2p+1 

implies easily 

n{z) = 

J=0 J=P+1 

It follows that 

^ _ r if .7=0,l,...,p, 

~ I A^-Pu) if j=p+l,p + 2, ...,2p. 

Therefore the above assumptions are equivalent to the following relations for the 
operator A and the generating element u, u € dom A^: 

(Ip) <=^ u e dom^P, 

(2p) ^ uedom^P, z{A{A - z)-^APu,APu) = 0{l), (4.1) 
(3p) uedom^P, z{A{A- z)-^APu,APu) + o{l) with a <^K.; 

in fact, in the last equivalence we have a — S2p+i. 
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Now we consider a function n e Nq with the property (Ip) for some p > 1. 
For < fc < p, by Sk wc denote the (A: + 1) x (fc + 1) Hankel matrix 

/sq Si • • • Sk ^ 

Sl 32 ■ ■ ■ Sfe+l 



Sk 



\sk Sk+l 



(4.2) 



S2k ) 



it is the Gram matrix associated with the fc + 1 functions to,ti, . . . ,tfc, and we 
introduce the Gram determinants 



Dk := det Sk = 



50 Sl 

51 32 



3k 
Sfc+1 



fc = 0,l,...p. 



(4.3) 



Sk Sfc+i • • • S2k 

Further, for k ~ 1, . . . ,p, TLk denotes the fc-dimensional subspace 

Hk ■■= span {to, ti, . . . ,tk-i} 

of 7i(n). Evidently, the subspace Hk is non-degenerated if and only if Dk-i ^ 0. 

In the rest of this section we suppose that Dp_i 7^ 0, that is, the subspace 
Tip is non-degenerated. If Dp ~ 0, then the function n with the given asymptotics 
(jp) is uniquely determined and rational of Mac Millan degree p, in fact, see [H 
pp. 22,23] 



n{z) 



dp{z) 
ep(z)' 



where the polynomials Cp of degree p and dp of degree p—l are defined below. To 
exclude this (simple) case we often suppose that even Dp 7^ 0; clearly, this implies 
Dp-i ^ 0. 

As a basis in Tip we choose a system of elements e 7i(n) = L'^{o'), k = 
0, l,...,j5— 1, which is obtained from the system to,ti,...,tp by the Gram- 
Schmidt orthonormalization procedure. This so-called system of orthogonal poly- 
nomials of the first kind, associated with the function n is defined by the following 
properties, j, fc = 0, 1, . . . ,p — 1: 

1. eo{z) EE 1/^s^, 

2. efe(z) is a real polynomial of degree fc with positive leading coefficient, 

3. (ej,efc) = 6jk- 
Then, see [fl (1.4)], 



ek{z) = 



^Dk-iDk 



So 

Sl 



Sfc-1 

1 



Sl 
32 

3k 
Z 



3k 
Sfc+1 



S2fc-1 
^fc 



, fc = l,2,...,p-l, (4.4) 
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and by this formula with k = p also a polynomial Cp can be defined. Evidently, 



Cp e Tip , and 



span {to, ti, . . . ,tfc} = span{eo, ei, . . . , efc}, fc = 0, 1, . . . ,p. 
The orthogonal polynomials ej, j = 0,1, ... ,p, satisfy the difference equa- 



tions 



bk-iBk-iiz) + afcefc(z) + 6^6^+1(2) = zckiz), fc 0, 1, . . . ,p - 1, 



(4.5) 



with real numbers Uk, fc = 0, 1, . . . ,p — 1, 6-1 = 0, and positive numbers bk, k = 
1, . . . ,p — 1, and the 'initial condition' eo(z) = 1/y^. Explicit formulas for ak, bk 
can be given, see [TI; we note that 



ao 



£1 

So' 



bo^ 



sq 



(4.6) 



The relation (|4.5p implies that with respect to the basis eo,ei, . . . ,ep_i of 
the space Tip the compression Ap of the operator A to Tip is given by the Jacobi 
matrix 



f ao bo 

bo ai bi 



\ 





V 



ap-2 bp^2 

bp-2 Up-i I 



(4.7) 



and that 



Acp-i — bp-2ep-2 + ap_iep„i + 6p_iep. 
The latter relation means for the orthogonal polynomials 

bp^2ep^2iz) + ap_iep_i(z) + &p_iep(z) = zep_i(z), 

therefore, the eigenvalues of Ap are the zeros of the polynomial Cp. For later use 
we write the last p — 1 difference equations (j4.5p explicitly in the form 



bici + 0262 + 6263 
^^262 + 0363 + 6364 

^ bp-2ep-2 + dp-iCp^i + bp-It 



~ zei 
= ze2 

= 263 



(4.8) 



this system of homogeneous equations for eo, ei, . . . , Cp determines the orthogonal 
polynomials uniquely if we add the initial conditions 



I \ / \ z — ao 

eo[z) = -—, ei(z) = - — — ; 
V^o OoJso 



(4.9) 



the second condition is just the first equation in (14. 5p . 
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The polynomials of the second kind, associated with the function n € Nq, are 
the functions dk, k = 0,1, . . . ,p, defined as foUows: 



dkiz) = y/s^ 



ek{z) - efc(-) 



-,eo = 



efc(z) - efe(-) 



k = 0,l,...,p. (4.10) 



Hence do{z) ~ and dk is a polynomial of degree fc — 1, fc > 1. The definition of 
dk and the relation (|4.5p imply that 

bk-idk-iiz) + akdk{z) + bkdk+i{z) = zdk{z), fc = l,...,p-l. (4.11) 

Therefore the polynomials Cfc and dk satisfy ioi k — 1,2, . . . ,p — 1 the same differ- 
ence equations but with different initial conditions: 

(4.12) 



do(z)^0, di(z) = ^. 



For later use we write the difference equations (|4.1ip in the form 

flidi + bid2 = zdi 

bidi + a2d2 + 62^3 = zd2 

62^2 + a3C?3 + bsdi = zds 



(4.13) 



k = 0,l,...,p-l. 



^ bp-2dp~2 + o-p-idp^i + bp-idp — zdp-i. 

For any two solutions uq, . . . ,Up and vq, . . . ,Vp of the difference equations 
dm with 6_i = 0: 

zuk{z) = bk-iUk-i{z) + akUkiz) + bkUk+i{z), 

(vkiC) = bk~iVk~i{0 +akVk{C) +bkVk+iiO, 

the Christoffel-Darboux formulas hold: 

p-i 

'^{z~C)uk{z)vk{0 = bp-i{up{z)vp^i{C) - Up^i{z)vp{C)) ^^ -^^^ 

k—7n 

-fe,„_i(u„(z)u„,_i(C) - u,„_i(z)w,„(C)); 

in particular, 

dp(z)ep_i(2:) - ep{z)dp-i{z) = . (4-15) 

Op-l 

2. In this subsection we assume that n G Ng satisfies the assumption (Ip) for 
some p > 1, and we consider its Schur transform n from (|2.7p . For the following 
lemma see [51 Lemma 2.1], we sketch the proof. 



Lemma 4.1. Suppose that n satisfies (Ip) for some p > 1 

, . So Sl S2p 

n{z) = 2 2^ + 

i/ien «is Schur transform n satisfies (Ip-i) : 

- , . ,Sn Si S9,,_9 



1 
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with for j = 0, 1, . . . , 2p — 2 



51 So 

52 Si 
Sj+2 Sj+i 





So 
Sj_i 








51 So 

52 Si 



(4.16) 



Proof. Write 



n{z) — z I 1 + 



Sl 



S2 



soz sqz^ 



S2p 
2p 



+ 



Sqz 



z + 



Sl 
So' 



If we set 



then 



S2 



S2p 



1 



1 + q{z) 



Sqz SqZ^ 



= 1 - qiz) + 



SqZ 



2p 



1 



2p-l 



2p 



1 + q{z) 



,2p-l 



r2p-l 



and 



n{z) 



l + g(z) 



si\ l + g(z) 



So/ 1 + 



Sl , 



-z,(.) + -(l + ,(z))^^^ 

'2 , , S2p / 1 



soz 



soz 



2p 



S2 

soz 



S2p 



saz 



2p-l 



z'^P 

1 



r2p-l 



f —q{z] 
So 

— 

So 



1 



l + qiz) 
1 

1+^' 



which is of the needed form. Formula (|4.16[) for the coefhcients Sj can be obtained 
by equating powers of z from both sides of the equahty 



n(z) [ z ^- + niz) 

So 



-So. 



□ 



Now we can formulate the main result of this subsection. 

Theorem 4.2. Let n e No satisfy condition (Ip) for some p > 2. If e^ and 

dfc, k = 0,1,..., p, denote the polynomials of first and second kind associated 
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with the function n, and Ck and dk, k = 0, 1, ... ,p — 1, denote the polynomi- 
als of first and second kind associated with the Schur transform n of n, then for 
A: = 0, 1, . . . ,p — 1 the following relations hold: 



ek{z) = -^dk+i{z), 



(4.17) 



dk{z) 



soek+i{z) + 



'So 



bo 



Ck+iiz) - Ck+li-) 



So 
ei ) = 



1 



dk+i{z) 

ek+i{z) - efe+i(-) 



ao 



(4.18) 



Remark 4.3. (i) Here we write e and d for the polynomials of first and second kind 
associated with the Schur transform n of n, but the reader is reminded that these 
functions are not the Schur transforms of the polynomials e and d. 

(ii) If efc and d^ stand for the polynomials of first and second kind associated with 
the Nevanlinna function — 1/n, then for fc = 0, l,...,p— 1 



efc(z) = (ife+i(z), (Ifc(z) = -efc+i(z) + 



1 

So 



£l 
So 



dk+i{z). 



The first equality readily follows from the fact that the spectral functions of n and 
— 1/n only differ by a factor sq. The second equality can be obtained by tracing 
the proof below; the only difference lies in (|4.19p : with evident notation, it should 
be replaced by 



So 



S0S2 - Sf 



Proof of Theorem H.'Il For the function n, again with evident notation, we have 



So 



S0S2 



(4.19) 



and, as a consequence of Theorem 13.1 



/So ^0 J3 

60 ai bi 



yo 








bp-z 



For the e^, j = 0, . . . ,p — 1, wc find 



\ 


bp-z 
ap-2) 



/ai 61 

61 02 62 



Vo 








\ 




ap-2 fep-2 
bp-1 Op-i/ 



(4.20) 
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and 61,62,. 



>ep-i 



follow from the equations (see (|4.5[) ): 



oieo + 6161 
6160 + 0261 + 6262 
^261 + 0362 + ^363 



= zeo 
= zei 

= 262 



Since these equation coincide with (|4.13|) we obtain 

ej = cdj+i, j = 0, 1, . . . ,p- 1. 

The constant c can be determined from the initial condition eo = c di , which gives 
6=1/ ^/so. Therefore 



1 



dj+iiz), j = 0, 1, . . . ,_p- 1, 



and (|4.17p is proved. 

For the polynomials of second kind dj we obtain in a similar way 



do(z) = 0, di{z) 



So 
bo ' 



aidi + 61^2 

bidi +02^2 + 62^3 

62^2 + a3c'3 + bsdA, 



zdi 
zd2 
zds 



Dp-3Mp-3 



a. 



p-2"p-2 + Op-2«p-l 



p-2 



(one equation less than in (j4.13p ). These equations can be written as 



02^1 + 626^2 

62^1 + 03^2 + b^d^ 

63^2 + 04^3 + b4d4 

^ bp^2dp-3 + ap^idp^2 + bp^idp-i 



zdi 
zd2 
zds 



zdp-2. 



(4.21) 



The last p — 3 equations of this system coincide with the last p — 3 equations of 
(j4.8p and (|4.13|1 . Therefore a solution vector {dj of the last p — 3 equations of 
(j4.2ip can be obtained as a linear combination of the solution vectors (6^)2 and 
{dj )2 of the last p — 3 equations in (|4.8p and (|4.13p : 



dj = 76^+1 + Sdj 



+ 1: 



j = l,2,...,p-l. 
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Now 7, 5 have to be found such that these relations hold also for j = with 



(io(z) = 0, and for j ~1 with di{z) ~ — — . Since ^0(2:) = it follows that 

^0 



= 7ei(z) + 5di[z) = 7 



which is satisfied for 



z — flo 1 



v/fo 



X i - si 
-£, = e \ 5- 



So 



60 



So 



So . 



The relation di{z) — " implies 
^0 



di{z) 



-62(2;) + ( — ] d2{z] 

\sq So J 

(z — ai)ei — 6060 z — ao {z — ai)di 

/ 2:-ao 1 

/ {z-ai)— 



^sq 



So 



bo 



bo ' 



flO Oq 



61 



So 



61 



hence 



&o 



^iv/so 



ysp 

So 



60 
foo^ 



h 

bo 



According to (|4.20p we find e = ^/sq- This proves the first equality in (|4.18p . The 
remaining equalities follow from (|4.10p and the second equality in (|4.9p . □ 

3. In this subsection we give a second proof of Theorem 14. 21 using asymptotic 
expansions, see [1] (1.34b)]: Assume that n E Nq satisfies (Ip) for some p > 2, 
that is, 



then 



n{z) = 
dp{z) 



So 

z 



£1 

r2 



S2p 



£0 

Z 



£1 

„2 



S2p-1 



+ 



1 



o 



1 



72p+l 



(4.22) 



ep(z) z z^'P 

According to [H the second to last formula on p. 22] the function —dp/cp is a 
Nevanlinna function and by |131 Bemerkung 1.11] there is a real number t2p such 
that —dpi Cp has the asymptotic expansion 

dp{z) _ £0 _ £1 S2p-1 

z -''^ 



ep{z) 



y2p 



t2p 



1 

^2p+l 



(4.23) 
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By Lemma [4. 1[ 
and hence by (|4.22p 

dp-ljz) _ So Si S2p-3 

By Lemma |4.1[ the Schm- transform of the function —dp/cp in (|4.23p has the 
asymptotic expansion 

r{z) 



^\-^]- (4-24) 



dp' 




_ soep(z) 


- 2 H 




dp{z) 


So 




fo 


Sl 


S2p-3 




z 


22 


^2p-2 




So 


Sl 


S2p-3 




z 


22 


^2^-2 



■ dp{z) 

^2p-2 

_22p-l ' " \ 22P^1 



(4.25) 



2p-l 



where only the number ^2^-2 depends on t2p according to formula (|4.16p . Here the 
polynomial r, defined via the second equality sign, is given by 

r{z) = soep(2) ^ (^'^ ^p(^) 
and its degree is < p. Comparing (|4.24p with (I4.25p . wc find that 

r{z) dp-i{z) I 



2p-l 



dp{z) ep-i{z) 
The degree of the product dpCp-i equals 2p — 2 and hence 

r{z) _ dp_i{z) 
dp{z) ep_i(2)' 
which readily implies that for some number fc 7^ 

ep_i(2) = kdp{z), dp^i{z) = k ^soep(2) ~ (^^ ^ ~^ dp{z)^ . (4.26) 

We claim k ~ 1/ -^/sq . With the proof of the claim the proof of the theorem is 
complete. 

To prove the claim we note that the leading coefficient of the polynomial 
is equal to y/Dk-i/Dk and that, by (|4.5p . 



Dk 1 Dk 



Dk+i bk V Dk 



Hence 



Dp-i _ 1 1 [Do 1 f f 1 



Dp 6p_i bi V Di &p_i bi 60 ^/s^' 
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and, similarly, because of (|4.20[) and (|4.19p . 



lDp_2 _ 1 1 1 _ 1 11 



Y Dp^i 6p_2 ^0 Vso ^p-i bi bo 
From (|4.26[) we obtain 



1 Dp^2 ,. dp{z) ep{z) Dp-i 

- \ / ^ — = lim = So lim — = So A / — — , 

Dp^i zP-^ ".-oo zP ]] Dp ' 

that is, 

1 11 11 11 



bp-i bibo kbp^i bibo 
Therefore, k = I/^/sq and the claim holds. 



5. Reduction via a p- dimensional subspace 

Let again n E Nq with the property (Ip) be given. We decompose the space 7i(ri) 
with Hp = span {to, ti, . . . ,tp_i} as follows: 



(5.1) 



Then, evidently, Cp G TL'p. The corresponding matrix representation of the operator 
A is 



(5.2) 



Ao B 
B D 

with Aq given by the Jacobi matrix ^o from (|4.7p . 

B ^ &p-i( ■ , ep_i)ep, B = 6p-i( • , ep)ep_i. 

The operator Aq is bounded and self-adjoint in Hp, and D is self-adjoint (possibly 
unbounded) in Tip. In the following theorem we express the function n by means 
of the entries of the matrix in (15.21). We set 



aoo{z) ■■= {{Aq - z) ^M, u), aii(z) := ((^0 - ^^) ^ep_i,ep_i), 

a{z) := {{Aq - z)~'^u,ep-i) = {{Aq - z)~^ep_i,u), (5.3) 



and 



f(z) 



000(2) a{z) 
a{z) aii(z) 

The last equality in (|5.3p follows from a{z*)* ~ a{z), in fact, by Cramer's rule, 
o-{z) — v^/ ('^o — 2) if p = 1 and a{z) = (— l)^~^i/so 60 . . . 6p_2/det (^0 — z) if 
p > 2. 
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Theorem 5.1. Suppose that the function n G Nq satisfies for some p > 1 one of 
the assumptions {jp),j — 1,2,3, and that Dp ^ 0, see (j4.3p . Then 

f{z)np{z) - apoiz) 

n{z) = — — , 5.4) 

aii(z)np(z) - 1 

where np{z) := {{D — z)~^Up,Up) with Up := bp-iCp. The function Up belongs to 
No and satisfies the assumption (jo). Moreover, for k > 1 we have 

u e dom A^^'' Up e doiJiD'^, (5.5) 

and n satisfies the assumption (jp+fc) if and only if Up satisfies the assumption 

Remark 5.2. If the operator representation (|3.2[) of n is n{z) = (^{A— z)~^u, u) with 
the space Ti ~ Ti-ln), then, according to Theorem l3.1i the operator representation 
of Jii = n is given by 

P-H' An 
\\u\\ 

the operator representation of n2 = ni by 

Pn'A^u 



M\\u,\V 

and, via induction, the operator representation of np = np_i by 

P-H' APu 



rL'2^ne{u,Au}, A2 = Pm^Aiw^, 

iction, the operator representation of 
H; = 7ieKA.,...,^-M, Ap = Pn,A\n>^, ^p HK|| . . . • 
Note that = ^fsa and, by Theorem 15.11 \\uj\\ = bj-i, j ^ 1,. . . ,p. 
Proof of Theo7^em \5.l\ With the matrix (|5.2p . the equation {A — z)x ~ u becomes 
{A(3 - z)xi + bp^i{x2,ep)ep-i = u, (5.6) 
bp-i{xi,ep-i)ep + {D - z)x2 = 0, 

where x is written as a; = [xi ^2)^ according to the decomposition (|5.ip . The 
second equation imphes 

X2 = -fop-i(a;i,ep_i)(D - zy^Cp. 

We insert this into (|5.6p and apply {A^ — z)^^ to get 

xi = (Ao - z)~^u + 6p_i(xi, ep„i) [{D - z)~^ep, Cp) {Aq - z)~'^ep-i. (5.7) 

Now take the inner product with ep_i and solve the obtained equation for (xi, ep_i): 

, ^ ^ ((Ao - z)-^M,ep^i) 

ixi, ep_ij ^ _ ^2^^ _ ^^^^ _ • 

Observing that n{z) = ((A — z)~^u, m) = (xi,m), the relation (|5.7p yields (|5.4p . 
To prove (|5.5p . denote by P' the orthogonal projection onto H.'^ in (|5.ip . If A: — 1, 
then u G domA^"*"^ is equivalent to v :— A^u G domA. Since D is the only 
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entry in the matrix of (j5.2[) which is possibly unbounded (B and B are even one- 
dimensional), V € A is equivalent to the fact that the non-zero component P'v, 
which is a multiple of Cp, belongs to dom_D. If fc = 2 we observe that 



( Al + BB A„B + BD\ 

' BAo + DB BB + I ' ^ ' ^ 



Since ep £ domZ? the operators DB and BD and hence all the entries in the 
matrix representation of A^ except possibly are bounded. Now u G domA^'+^ 
is equivalent to v = A^u e domA^, and hence, by (|5.8p . P'v E douiD^. The claim 
for arbitrary k follows by induction. 

The last claim of the theorem for j ~ 1 follows immediately from (jS.Sp and 
the first equivalence in (|4.1|) . For j = 2, 3 we also use the equivalences in (|4.ip . A 
simple calculation yields 

with 

Ruiz) = Si{z)-\ i?i2(z) = -bp^i {{D - z)-i • , ep) Si{z)-^ep.i, 

R22{z) = {D- z)-i + bl_, (5i(z)-iep_i, ep_i) {{D - z)-^ ■ , Cp) (D - z)-^ep, 

where Si{z) := — z — [{D — z)~^ep, ep) ( • , ep_i)ep_i, the first Schur com- 
plement. It is easy to see that for f,g€ Tip we have 

lim iy{Si{iyr^f,g) ^ ~{f,g). 
Now we observe the relation 



A{A-z)-^ 



%i?ii(z) + BRi2{z*)* Aoi?12(z) + BR22{z) 

^BRn{z) + DRi2{z*r BRi2{z) + DR22{z) 



and the fact that for z = iy, y —* oo, for example zAoRii{z) = zA()Si{z) has a 
limit, 

zAoRi2{z) = -zbp^i{{D-z)-'^-,ep)AaSi{z)-^ep-i=o{l), 

zDR22{z) = zD{D~z)-^ + zhl_^{Si{zy^ep,ep){{D-z)-^-,ep)D{D~z)-~^ep 

= zD{D - z)-^ +o{l), 

etc. These relations imply for example with v — AP'^^u 

z {A{A - z)-^v, v) ^z {D{D - z)~^P'v, P'v) + 0(1). 

Since P'v £ span{ep, Cp+i, . . . , Cp+fc}, ep, Cp+i, . . . , ep+fc_i e domZ? and hence 

z {D{D - zY^x' ,x') = 0(1) for x' e span{ep, Cp+i, . . . ,ep+fc}, 

and since P'v has a non-zero component in the direction of ep^k the claim follows 
from (SHI). □ 
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6. Representation of the transformation matrix by orthogonal 
polynomials 

The 2x2 matrix function which generates the fractional hncar transformation 
(|5.4[) we denote in the foUowing by V: 



V{z) 



1 / f{z) -000(2:)^ 
aiz) \aii{z) -1 , 



(6.1) 



In this section we express V by the polynomials of first and second kind. To 
this end, the elements of Tip are considered as column vectors with respect to the 
basis Co, ei, . . . , Cp-i. 

First we solve the equation (^0 — z)x = Cp^i in Tip. With the Jacobi matrix 
^0 from (|4.7p this equation becomes 



or 



{Ao - z) 



\ J V 1 y 

According to the definition of the orthogonal polynomials of first kind, the solution 
of the system with the 1 in the last component of the vector on the right hand side 
replaced by —bp-iep{z) is the vector with components 60(2:), ei{z), . . . , ep_i(z). It 
follows that 

ej„i(z) 



fep_iep(z) 



, j 1,2, . . . ,p, 



and hence 



{{Aq-z) iep_i,ep_i) = (a;,ep_i) = - /^^ 

bp^iep(z) 



that is. 



{{A-z) ^ep_i,eo) = (x,eo) 



011(2) = -7 TT' ^{^) = -- 



£0(2) 
5p_iep(2)' 



1 



, , y -^v-. , , V (6.2) 

&p_iep(2) &p_iep(2) 

Next we solve the equation (^0 ^ z)x = u ~ ysoeo. As above, in matrix form it 
becomes 



{Ao - z) 



6 _ 

Up / V ; 
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According to the definition of tlie polynomials of the second kind and because of 
bo di{z) ~ we have 



/ \ 

di{z) 



(Ao - z) 



( 



dp-2{z) 

V dp-i{z) ) 
It follows that 











6 



and hence 



= (A - z)-^ 



\ -bp-idp{z) J 



f V^o\ ( 














( \ 
























- bp^idp{z) 




















v 





) 




1 1 / 



V ; 





di{z) 



dp{z) 
ep{z) 



V dp-iiz) J 



( eo(z) \ 
ei(z) 

\ ep-i(z) / 



am{z) ^ {{Ao ~ z) ^u,u) = - 



dpjz) 
ep{z) ' 



(6.3) 



Inserting the expressions from (j6.2p and (j6.3p into (|5.4p and observing the relation 
(j4.15p we find that V can be written as 



V(z) = f^'^P-i(^) -bp.idp{z)\ 



and hence we obtain the following theorem. 

Theorem 6.1. //, for some integer p > 1, the Nevanlinna function n satisfies one 
of the assumptions {jp), j € {1, 2, 3}, and Dp ^ theM the following relation holds: 



n{z) = {{A z)-\, u) ^ yP-ii-)M^) + yidp{^) 

ep_i[z)np{z) + bp-iep{z) 

where np[z) = ((Z? — z)^^Up^Up) , Up = bp^iCp. 

Remark 6.2. (i) Using and we obtain from with p = 1: 

n{z) 



(6.4) 



So 



z — — + ni{z) 
So 

and hence, because ni{z) = o(l), ni is the Schur transform of n: ni = n. For p>2 
we obtain from (|6.4p and (|6.4p with p replaced by p — 1 and with the help of 
and (|^7TS)) that 



np_i(z) = - 



bl-2 



z — ap-i + np{z) ■ 



hence 



np_i(z) = - 



^P-ibp- 
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and rip = rip-i. 



(ii) From (EH), (|iT71) . ((^T^ . and g^Ol), we obtain 




{z - ao) 



dp^2{z)np{z) + bp-2dp-i{z) 
ep^2{z)np{z) + bp-2ep^i{z) 



Since, according to Theorem 14. 2 [ dk and Ck arc the polynomials of first and second 
kind associated with n , this formula implies that the function Up is the p — 1-st 
Schur transform of n. 

The 2x2 matrix polynomial V, which generates the fractional linear trans- 
formation (j6.4|) . has the property 



Therefore V{z) ^ exists for all z G C and we can form the polynomial matrix 



detV{z) = 6p_i {dp{z)ep-i{z) - ep(z)dp_i(z)) = 1. 



With 




V is J-unitary on the real line, that is. 



V{z)JV{z)* = J zeR. 



function 




with 




6p_i {dp{z) dp_i{0) 
bp-i {dp{z) ep^i{0) 
6p_i (ep_i(z)dp(0) 
6p_i (ep_i(2)ep(0) 



dp-i(z)dp(0)) , 
dp-i{z) ep{0)) , 
ep(z)dp_i(0)) , 
ep(z)ep_i(0)) . 



A straightforward calculation leads to the relation 



n 



(z) = ({A - z)~^u,w) = 



p['>{z)hp{z)+p\,'>{z) 
q[P\z)hp{z) + ql''\zy 



where 
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With relation (|4.14p . the following formulas can be obtained, compare [U 1.2.4]: 

p-i 

fc=0 

p-i 

(P)^z) = l + z^dfe(z)efe(0), 



Pi 



k=0 

p-l 



p-i 



fc=0 



7. Transformation by means of a w-resolvent matrix 

Given again a function n G Nq with one of the properties (jp), j ~ 1, 2, 3. Besides 
the decomposition (|5.ip we consider the decomposition 

7Y(n) = Hp+i e 7^", = span{7^p, e^} = span{eo, ei, . . . , Cp}, 

and in the space 7ip+i the restriction 

5 := A\n, = 



B 

This restriction is a non-densely defined symmetric operator in 7ip+i with defect 
index (1, 1), and. evidently, the given function n = (^{A — z)~^u, uj is one of the 
M-resolvents of this operator S. Hence n can be represented as a fractional linear 
transformation of some function g e No by means of the w-resolvent matrix W = 

W2l[z)g{z) + W22[Z) 

Such a u-resolvent matrix W can easily be calculated. To this end we fix 
a self- adjoint extension of S in 7ip+i, which means that we fix some 7 €E M in 
the right lower corner of the matrix representation of S with respect to the basis 
Co, ei, . . . , Cp of 7ip+i. Denote this matrix or self-adjoint extension of S in 7ip+i 
by Aq,^: 

According to [14j this w-resolvent matrix W is given by the formula 

^ 1 [ iR2u,u) {R2u,u)Q{z)~ {u,ip{z*)){ip{z),u) \ ,^2) 
^ ' iu,ip{z*)) \ 1 Qiz) y ' 
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where Kl = {Aq^^ — z)""^ , (/^(z) is a defect function of S corresponding to the self- 
adjoint extension Ag^-y, and Q is the corresponding Q-function. An easy calculation 
yields 



A{z) ' 



where 



A(z) = 2 - 7 + 6p_iaii(z). 
Since S = Ao^^lfi^ and hence, in terms of linear relations, 

S* ~ {{x, Aq^jX + Acp} I X G Hp+i, A e C}, 
it is easy to check that for (p{z) with {v(z), z(p{z)} G S* we can choose 

-1 



ip{z) = (Ao.j - z) 



1 



and then the Q-ixmction, which is the solution (up to a real additive constant) of 
the equation 

Q{z)-Q{(r 



becomes 



z- C* 

Q{z) = 



{^{z),v{0), 



A{zy 



Inserting these expressions into W from (|7.2p we find 



Wiz) 



A(z) 
6p_ia(z) 



aoo(z) ■ 



a{zf 
-A(^ 



^00 



V 



A(z) 
-1 
A(I) / 



(7.3) 



Observe that W{z) is J-unitary on the real line. Next we establish the connection 
between the matrix functions V from (|6.ip and W from (|7.2|) . in fact we find a 
simple expression for V~^W . We have 



V{z)-^ 



1 / -1 aoo(z) \ 

a(2) V-aii(z) aoo(2:)aii(z) -0(2)7 ' 
Multiplying this matrix from the right by W(z) from (|7.3p we obtain 





F(z)-il^(z) 



z — 7 



(7.4) 



Theorem 7.1. // the function n E Nq has one of the properties (jp), j G {1,2,3}, 
then the matrix functions V from (16. ip and W from (j7.2p are connected by the 
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relation (j7.4[) . Therefore for the Nevanlinna functions Up in (|5.4p and g in (|7.ip 
f/ie following relation holds: 

np{z) = (7.5) 

If l/g{z) = 0(1), then formula (|7.5p implies that rip admits the asymptotic 
expansion 

l^l^i , / 1 



np(z) ^ 

and — 1/5 is its Schur transform: ~l/g = rip. Hence the number 7, which defines 
the self-adjoint extension of S, corresponds to the number Op. 

Remark 7.2. If instead of a self-adjoint operator extension Aq.^ of S we choose the 
(multi-valued) self-adjoint relation extension of S: 

^0,00 = 5 + span {0, Cp] ^ Aq + span {0, Cp}, 

then wc obtain 

i?r = (^0,00 - z)-^ = (^0 - z)-^p, 

where P is the orthogonal projection in 7ip+i onto Tip, 

<<5(^) = ( I, Q[z)^ z + bp_^aii{z), 

so that 

W^(z) 



z ~ 








z 











1 




z 




1 


bp-i ^ 



bl-Mz) y 1 

V{z)-^W{z) = 
and instead of (j7.5p wc have 

2 + 5(2) 

Thus if g(z) = o(l), then Up has the asymptotic expansion 

6Li / 1 



np(z) = 

and g is the Schur transform of Up: g = rip. 



y2 



A more explicit form of the resolvent matrix W from (j7.2p can be obtained 
following [H] and [2]. To this end we decompose the space 7ip+i as 

T^p+i ~ ran(S' ~ z) + spanu, z e C, a(z) ^ 0, 
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(+ stands for direct sum) and denote for y e Ti-p+i by P{z)y the coefficient of u 
in the corresponding decomposition of y: 



y = {S - z)x+{P{z)y)u 



(7.6) 



with some x G doni5 — Tip. Further, define Q(2;) y = {{S — z) ^(y — {P {z)y)u) , uj 
Then, according to [13] , the resolvent matrix can be chosen to be 



W°{z) ^h + z 



Q(z) 
-P{z) 



(Q(0)* -P(0)*)J, J 



f 

-f 



(7.7) 



We derive an exphcit expression for W^{z) := W{z)W{Q)~^, following [5]. To 
this end, for the vectors and operators we use matrix representations with respect 
to the basis to( = u), ti, 

this basis. We denote by 6 the (p + f ) x (p + l)-matrix 



tp. Recall that Sp is the Gram matrix associated with 



6 



f 




Vo 











1 0/ 



and by e the first column in the {p + f ) x (p + f ) identity matrix. Then S and u 
correspond to 6 |cp+{o} ^i^d e. First we apply the operator (/ — z&*)^^ to (|7.6p 
and observe the relation 

e*(/-2e*)-i(e-z)x = 0, a;eCP + {0}. 

It follows that 

e*(/ - ze*)-^y = (P(z)2/) e*{I ~ z6*)-\ = P{z)y. (7.8) 
Further, observing that e* Sp = (sq si ■ ■ ■ Sp) we obtain 

Q{z)y = ((5-z)-i(2;-(P(z)2/)7.),ii) 

= e*5p ((S - z)-'y - (6 - 2)-ie (e*(/ - ;^S*)-iy)) 

-(so si ••• Sp)((e-z)-i(/-ze*)-(e-z)-iee*)(/-z6*)"iy 

= (so si ••• sp)e*ii-ze*)-^y 

= (0 So si ••• Sp_i) (/ - z6*)"^y, 
where for the second last equality sign we have used that 

(6 - z)-\l - z6*) -{6- zY^ee* = 6*. 
Together with dLS]) we find 



Q(z) 

-p(z; 



So si 
-10 
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and (17.71) becomes 



so 




Sp-l 



So 
Vp-i / 



1 

-1 



8. Reproducing kernel spaces: reduction via resolvent invariant 
subspaces 

In this section wc start from the operator representation of the Nevanhnna function 
n in the corresponding reproducing kernel space C{n) with kernel 

^ / X nlz) — n(w)* ^ ^ 1 ^ 

z — w 

see Section [21 (3). The operator A is introduced via its resolvent {A — z)^^ which 
is the difference-quotient operator Rz defined by 



(8.1) 



If n satisfies one of the assumptions (jp), then, by ^ Lemma 5.1], the fimctions 

./o(C) = n{C), /i(C) = C«(C) + so, . . . , /p(C) = CniC) + C-'so + ■■■ + Sp_i 

all belong to C{n) and 

{fk,fj)c{n)^Sj+k, j, fc = 0, 1, . . . ,p. (8.2) 

In particular, u := n Cz C{n), and by the reproducing property of the kernel L„ we 
have 

n{z) = {{A - z)~'^u,u)c{n}- 

By Uj we denote the class of all 2 x 2 matrix polynomials O which arc J- 
unitary on M and for which the kernel 

j-e{z)jQ{wr 

Ke(z,w} = 

z — w* 

is non-negative. The reproducing kernel Hilbert space with this kernel will be 
denoted by H(6); its elements are 2-vector functions. The matrix polynomials V 
and W considered in the previous sections belong to Uj: this follows from the 
Christoffel-Darboux formulas (|4?T4l) for V and from (fTi]) for W. Note that if 6 
belongs to Wj, then det 0(z) = c, where is c is a unimodular complex number, 
because the determinant detQ{z) is a non- vanishing polynomial in z. 

The following theorem was proved in [U Theorem 8.1], even in an indefinite 
setting. 
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Theorem 8.1. Let n G Nq and suppose that there exists a matrix polynomial 



such that the mapping 

u > (l — n) u 

is an isometry from 7i(0) into C{n). Define the function fi by 

^ a{z)h{z) + b{z) 
^ ' c{z)fiiz) + d{z)' 

Then the following statements hold. 

(i) fi is Nevanlinna function. 

(ii) The mapping g t-^ f: 

/(C) = (a(C)-n(C)c(C))5(C) 

is an isometry from C{n) into C{n). 

(iii) We have 

= (1 -n)n{Q)®{a-nc)C{n) 

and the mapping 



W : C{n) 9 / ( ^ ) e 



where /, u, and g are connected by the relation 

/(C) = (1 -"(0) u(C) + (a(C) - n(C)c(C))g(C), 

is a unitary mapping from C{n) onto TL{Q) C{fi). 
(iv) The mapping WR^W* is of the form 



WRM* 



Ruiz) Ri2{z)\ . fvie)\ fr{e)\ 

R2i{z) R22{z)) ■ \C{n) ) \C{n) ) 
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with 



Ruiz) =i?.-^(i?.e)(.)(^^^'^)(o l)E^ 

= -{a{z) - n{z)c{z))Ke{ ■ , z*) (Jj E,, 

R2i{z) =-;fc^(^-")(-)(o I) E, 

^ Lf.i-,z*){0 1) E,, 



k{z 

R22{z) =R,-^(R,n){-)E, 
k{z) 

where Rz is the difference- quotient operator, E^ is the operator of evaluation 
at the point z on any reproducing kernel space, and 

, / N / N w \ 7/ \ det Q(z) 

k{z) = c{z)fi{z) + d{z) - ^ ' 



a{z) — n(z')c[z) 

We mention that formula (|8.3p corresponds to the relation (|5.9p above. 

A space of functions is called resolvent-invariant if it is invariant under the 
difference-quotient operator R^ as defined in (|8.ip . In the following lemma, with 
a resolvent-invariant non-degenerate invariant subspace of a certain inner product 
space a 2 X 2 matrix function is associated. 

Lemma 8.2. Let Ai be a finite dimensional resolvent-invariant space of 2~vector 
polynomials endowed with an inner product ( • , • ) such that 

{RJ, g) - (/, R.^g) ~ {z - w*){RJ, R^g) = giwfJf^z), (8.4) 

and let that Mi he a resolvent-invariant non- degenerate subspace of Ai. Then 
there exists a 0i G hij such that 

(i) Mi = n{Qi), 

(ii) M = H(ei) © QiAf where Af ^ Q^^Mi is a resolvent-invariant space of 
2-vector polynomials, for which the relation ()8.4p holds if equipped with the 
inner product 

{Qi'f,Qi'g)M = {f,g), f,geMi^. 
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Relation (|8.4[) is often called de Branges identity, see [7] and, for ftirther 
references, [TT]. That J\f consists of 2~vcctor polynomials is due to fact that 
Q^^{z) = —JQi{z*)* J is a matrix polynomial. The other claims of the lemma 
follow from [SI Theorem 3.1]. 

Now we formulate and prove Theorem l6. ll again in the context of reproducing 
kernel spaces. 

Theorem 8.3. //, for some integer p > 1, the Nevanlinna function n satisfies one 
of ttie assumptions (jp), j € {1,2,3}, and Dp ^ then the following relation holds: 

^^^^^yp.^iz)npiz) + bp.,dp{z) ^^^^^ 
ep^i[z)np{z) + bp-iep{z) 

where fip is a Nevanlinna function such that fip G C{hp), and 

fip{z) = {R^fip,fip)c{np)- (8.6) 

Comparing (|6.4p and (|8.5p we find that fip{z) = np{z), the p-th element in 
the sequence obtained by applying the Schur transformation p times starting with 



Proof of Theorem 18.31 Let A4 be the linear space spanned by the p + 1 2- vector 
functions 

fo(C) = ,f:(C) - (!"^) ,...,f.(C) - (^°^'" V^,^''-') (8.7) 

and equipped with the inner product which makes the map u i— > (l — n) u an 
isometry from M into C{n), see (|8.2p . Note that if n has the integral representation 
the elements of Ai arc of the form 

irjR(fmdait) 

-/(C) 

where / is a polynomial of degree < p. Indeed, it suffices to show this for the basis 
elements of A^: If /(C) = C'' , then 

Refit) = ^i-^ = c^-i + + . . . + t^-^c + 



and hence, on account of (|1.3p . 

{Rcfm da{t) = soC^-i + + . . . + Sp_2C + sp- 

It follows that A4 is also spanned by the polynomial vectors 

j = 0, 1, . . .,p. 



-d, 



where Cj and dj arc the polynomials of first and second kind associated with n, 
see dOl) and (|jTTU)) . 
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Let M.p be the space spanned by the first p of the 2-vcctor functions in (|8.7p . 
Since Sp-i from (|4.2[) is a positive matrix, the space A^p is non-degenerate and 



M = Mp e span 



-dp 



As both Ai and Aip are resolvent-invariant spaces, by Lemma [8.21 we have that 
for some 9i G Uj, which is normahzed by 6i(0) = I2 (and hence det6(z) = 1), 

Mp = H{ei), M=n{Qi)®QiN. 

Here A/" is a one-dimensional resolvent-invariant space, which, when equipped 
with the induced inner product, satisfies the de Branges identity and therefore 
is spanned by a constant J-neutral vector {a /?) such that 



0i(^) ( ) = bp-i 



'dp{z) 

For A G R denote by C\ the constant J-unitary matrix 



\a a 



a = 0. 



A/3 

Then there exists a A such that 

^a{z) -bp^idp{z) 
^c(z) 6p„iep(z) 

where a and c are polynomials such that dcg c < deg Bp — p. The inclusion 



e(z) := Qi{z)Cx = 



Ro en{e) = n{ei) = Mp 

implies that deg a < p — 1. From det 0{z) = 1 it follows that x = a and y = c are 
polynomial solutions of the equation 

x{z)ep{z) + y{z)dp{z) = 

Op-l 

Since all polynomial solutions of this equation are given by 

x{z) ~ a(z) — s{z)dp{z), y{z) — c{z) -f s{z)ep{z) 

with some polynomial s, the solutions x = a and y = c have minimal degrees and 
because of that they are unique. Observing (|4.15p . we find 

a{z) = -dp-i{z), c{z) = ep-i{z). 

Hence 

6(2) = f-^P-i/^) -bp_idpiz)\ ^ ^ 
\ ep_i(z) bp^iep[z) J 

and Cx — 0(0) is the coefficient matrix of the fractional linear transformation 

(IS3D. 
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Define the function fip by (|8.5p . Then, according to Theorem l8.ll it is a Nevan- 
linna function. We show that fip G C{np). The function /p(C) = (l ^'^(0) fp(C) 
belongs to C{n) and, according to Theorem 18 . 1 1 fiii) . it can be written as 

(1 -n(C)) ipiO = (1 -n{0) u,(C) + (a(C) - n(C)c(C))gp(C) 

with Up £ Mp, gp € C{np), and the two summands on the righthand side are 
orthogonal. This orthogonality and the isometry of the mapping (l —nj imply 
that (0 7^) fp — Up e Mp and hence there is a non-zero complex number 7 such 
that 

Therefore 

(a(C) - n(C)c(C))5p(C) = -7 {dpiO + epiCMO) 

and 

-71(C) c(C)+a(C) 

Hence hp € C{hp). Equality (|8.6p follows from item (3) in Section [3] □ 
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